METHOD OF AND APPARATUS FOR THE ELIMINATION OF THE EFFECTS OF 
INERTIAL INTERFERENCE IN FORCE MEASUREMENT SYSTEMS, INCLUDING 
TOUCH-INPUT COMPUTER AND RELATED DISPLAYS EMPLOYING TOUCH FORCE 
LOCATION MEASUREMENT TECHNIQUES 

The present invention relates to the use of force 
measurement techniques, as in touch screen or other touch 
input computer cathode ray tube (CRT) or related displays 
to locate the points of touch on the face of the display, 
and under adverse circumstances where the display appara- 
tus is rather heavy and imperfectly or poorly or 
insecurely supported or mounted so as to be subject to 
inertial interference effects, such as vibration or sway, 
that may cause errors in the force measurement and result- 
ing computed touch location. 

While directed, accordingly, to obviating such iner- 
tial interference effects and rendering the force measure- 
ments immune to the disturbing effects of movement of the 
mass supported by the force measuring apparatus, as from 
such unsteady mountings or supports for the force measur- 
ing apparatus, the invention even more broadly finds use 
also with analagous problems arising in other force 
measuring applications, as well; such as, for example, 
those of weighing objects on a relatively insecure or 
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moving support, or effecting such weighing in Che presence 
of a relatively large tare mass or counterweight that is 
excited by rapid operation. , 

From still a broader viewpoint, the invention relates 
to enabling the accurate measurement of rapidly changing 
forces in situations where the associated mass displace- 
ments can cause significant inertial interference. 

BACKGROUND OF INVENTION 

As above noted, the invention has important applica- 
tion in touch input computer and 'related displays that 
employ touch force location measurement techniques. In 
particular, it provides significant immunity to inertial 
interference effects in specific systems such as those 
described in my earlier copending U.S. patent application 
Serial No. 857,241, filed January 25, 1992, for METHOD OF 
AND APPARATUS FOR TOUCH - INPUT COMPUTER AND RELATED DIS- 
PLAY EMPLOYING TOUCH FORCE LOCATION EXTERNAL TO THE DIS- 
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Before proceeding, however, to address such specific 
applications, it is believed helpful to consider first the 
general problems of the vibration, sway or other movement 
of heavy or poorly or insecurely mounted force measuring 
systems that may introduce errors into the force measure- 
ments and the resulting computations stemming therefrom* 

There is always some amount of mass associated with 
the "measurement side" of a force measuring device to 
which an unknown force is applied, including all or some 
portion of the measuring device itself, and the attachment 
and any related objects between them. If and when the 
device support changes motion, this supported mass is 
generally also accelerated. Since the force creating this 
acceleration passes through the device, the device output 
can no longer accurately reflect the unknown force. 

Additionally, the attachment, the device itself, and 
the device support, deform in some degree when the unknown 
force is applied; usually, in an elastic manner. Thus the 
supported mass must undergo some acceleration as the 
unknown force changes, at which times the device output 
responds not to the unknown force alone, but also to the 
inertial reaction force of the supported mass. 
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These effects are of practical concern when, for 
example, a force must be measured in the presence of' 
an unsteady support and an appreciable supported mass, or 
a rapidly changing force must be measured in the presence 
of appreciable supported mass and compliance. 

An approach to try to avoid such difficulties may 
reside in attempting to avoid the combinations of condi- 
tions listed above. Where such combinations are, however, 
unavoidable, two prior art techniques for amelioration are 
noted: viscous damping of the mechanical system itself, 
and linear filtering of a measurement signal derived 
therefrom. Damping is primarily applicable to situations 
such as weighing, where a force, once applied, remains 
constant until its measurement is completed. During such 
time as the force remains constant (exclusive of the 
inertial effects)., the intertial error decays exponent- 
ially. When the force must be measured as it changes, on 
the other hand, linear filtering may be particularly 
applicable. The filter may be of analog electronic, digi- 
tal electronic, or even mechanical construction, and is 
typically of low-pass or mixed low-pass and notch design. 
It passes low frequencies unaltered, but blocks those at 
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and/or near the resonant frequencies, of the suspended 
mass. A measurement is thus obtained which accurately 
reflects the low frequency part of the unknown force, but 
at the expense of abandoning measurement of the high fre- 
quency part. Both damping and filtering typically provide 
faster measurement than can be acheived with no ameliorat- 
ing technique: but not so fast as could be obtained were 
there no inertial interference present, or if there were 
some way to estimate such interference sufficiently 
ac cur a t e ly so that it could be subtracted from the total 
signal. 

It is to the provision of a novel ameliorating 
approach that the invention is directed, and will be here- 
inafter described in terms of the illustrative system of 
said copending patent application, wherein force measure- 
ment is used to locate finger touches on a typical CRT 
display monitor. Such application necessitates the accu- 
rate measurement of as much as 80 pounds, and must contend 
with flexure in the monitor housing and stand, and the 
table beneath, as well as in the force measuring platform 
itself. It has been found that good results (location 
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errors less than 0.1 inch) can routinely be achieved with 
such a system, provided the table is solid and the monitor 
housing rests on the platform through small, corner- 
located feet. The lowest resonant frequencies of such 
combinations run in the range of 6 to 10 Hz, and a linear 
filter that strongly blocks frequencies above 5 Hz has 
been found to provide adequate speed of response. 

Such is not, however, the case for less secure moni- 
tor supports. Many monitors, indeed, are supplied with a 
"tilt-swivel" stand. This form of support raises the 
monitor center of gravity and softens its support, parti- 
cularly with respect to fore-aft and side-to-side nodding 
motions. These two changes directly affect the two lowest 
frequency eigenmodes, dropping them to about 
3Hz for a large monitor on a tilt-swivel stand. A com- 
parable depression* of lowest resonant frequencies is seen 
when a free-standing folding table is used below the plat- 
form and monitor. In addition, disturbances such as bump- 
ing or leaning upon an insecure table can introduce large 
low-frequency inertial errors independent of the free 
resonant frequencies of the system. 
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While linear filtering with a lower cutoff frequency 
could, in some instances, be used to eliminate excessive 
interference from these lower frequencies, this would be 
only^at the cost of a touch response delay which is 
unacceptable. On the other hand, with only the higher 
frequency filtering in place, position errors can grow to 
as much as several inches for touches applied when the 
table and/or monitor are shaking from moderate bumping, or 
from several rapidly delivered strong touches. 

As before stated, it should be noted that analogous 
difficulties may arise in other force measuring 
applications, such as that of weighing objects on an 
insecure or moving support, or weighing in the presence of 
a large tare mass which is excited by rapid operation. 

The invention therefore provides a new solution to 
this class of problems. 

OBJECTS OF INVENTION 



It is accordingly an object of the invention, 
therefore, to provide a new and improved method of and 
apparatus for measurement of forces that are immune to 
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disturbing effects of an unsteady mounting or support for 
the force-measuring apparatus through canceling or elimin- 
ation of errors caused by such effects in the force 
measurements . 

It is a further obiect of the invention to provide 
such a novel method of and apparatus for measurement of 
forces that are immune to the disturbing effects of move- 
ment, in response to changes in the applied force, of the 
mass supported by the force-measuring apparatus, as well. 

Another object is to provide these immunities in a 
manner that can be conveniently applied to a wide range of 
circumstances, with characteristics of support and of 
supported mass which are not beforehand known, but for 
which the force-measuring apparatus can be easily field- 
calibrated by the user. 

Still another object is to provide an improved 
touch-input computer or related display employing touch 
force location, including externally to the display, that 
is rendered immune to inertial interference effects. 

Other and further objects will be explained herein- 
after and are more particularly delineated in the appended 
claims . 
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SUMMARY OF INVENTION 

In summary, from one of its broader aspects, the 
invention embraces in a touch-input computer and related 
supported display employing touch force location measure- 
ments, a method of eliminating the errors that may be 
introduced into force and/or torque measurements by 
undesired inertial interference motions of one or more of 
the support, mechanical system of the display and/or force 
measuring apparatus itself, that comprises, sensing one or 
more components of force and/or torque applied to the dis- 
play by touch forces to provide force and/or torque 
measurements uncorrected for inertial interference motion 
effects that may arise; sensing lineal and/or rotational 
acceleration of the display in response to such inertial 
interference motions; and correcting the uncorrected force 
and/or torque measurements in response to the acceleration 
sensing to reflect ^substantial elimination from the 
measurements of the effects of such inertial interference. 

Other features of the invention and preferred and 
best mode designs are hereinafter set forth. 
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Underlying the invention, however, are three related 
aspects, corresponding to the above-stated objects. 

1) Particular accelerometers are incorporated into 
the force-measuring device in sufficient number and manner 
to record all relevant degrees of freedom of motion of the 
support (if attached to the support side), or of the 
attachment to the unknown force (if attached to the 
measurement side). A typically different Linear combina- 
tion of these acceler ome t er channels is added in turn to 
each channel of uncorrected force measurement, with co- 
efficients of combination chosen such that the resulting 
sums reflect the desired force measurements, free of 
error, in the low frequency limit, or in the limit of an 
arbitrarily stiff connection to all of the supported mass, 
due to tilting or change of motion of the support. 

2) In addition, or alternatively if desired, sets of 
correction channels are derived corresponding to time 
derivatives of the uncorrected force measurements, includ- 
ing typically at least the set of the second derivatives, 
but possibly also sets of one or more derivatives of other 
orders. A typically different linear combination chosen 
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such that the resulting sums reflect the desired force 
measurements, free of errror due to changing the supported 
mass, to the extent that the support itself does not 
change motion. 

3) Coefficients of combination for correction 
channels are determined empirically, in situ. Such 
correction channels may comprise accelerometer readings, 
and/or time derivatives of various order of these or of 
the uncorrected force readings, and/or related quantities 
of similar value, such as quantities effectively approxi- 
mating the preceding over a desired frequency range, or 
linear combinations of any measured data containing the 
desired acceleration and/or derivat ive data* The force™ 
measuring device is installed into its normal operating 
context, which may vary considerably from installation to 
installation. To calibrate, data are then gathered in the 
same manner as for normal force measurement purposes, but 
while the applied unknown force is allowed to remain at 
zero, and then while the system is intentionally disturbed 
in various ways. The idea is to gather a reasonably rich 
collection of measurement sets from moments when forces 



-12- 



caused by acceleration of the supported mass are the only 
ones applied to the measurement side of the device. Dis- 
turbances may consist of pushing or shaking the system 
support in differing directions and patterns while data 
are gathered and/or pushing or shaking the supported mass 
at different points and in differing directions; data 
being gathered after the disturbing force is removed. A 
correction matrix may now be found, by methods later 
described, which relates the data collected from the 
correction channels to that gathered from the uncorrected 
force channels. The elements of this matrix comprise the 
desired coefficients of combination to achieve the above- 
stated corrections. 

The invention contemplates that any of these three 
independently novel aspects may be separately used, or 
that they may be used in various combinations. The parti- 
cular preferred application and embodiment described here- 
in benefits indeed from incorporating all three. In this 
manner, the method and apparatus of my said copending 
application is significantly improved to work accurately 
and reliably with a wide range of massive and softly 
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supported computer display monitors. A simple one-time 
calibration procedure is performed at display installation 
time, requiring no equipment, and being effectable in less 
than one minute to complete. For installations that would 
otherwise be plagued by large inertial effects, the typi- 
cal performance improvement achieved is a factor of about 
thirty times reduction in error of reported touch loca- 
tion. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will now be described with reference to 
the accompanying drawings, 

Fig. 1 of which is a side elevation showing an 

illustrative application in which the preferred 
embodiment supports a computer CRT display monitor 
for purposes. of user touch location as described in 
my said copending patent application; 

Figs. 2a and 2b illustrate the design of an appropriate 

accelerometer for use with the invention, in plan and 
exploded views, respectively; 
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Fig. 3 is an isometric view showing the disposition of 
components within the force-sensing platform, 
employing the accelerometer improvement; 

Fig. 4 is a simplified schematic diagram of electronic 

conversion and calculating circuitry suitable for the 
system of the invention; 

Fig. 5 provides a flow-chart of the firmware operation of 
the preferred embodiment. 

Fig. 6 illustrates the definition of vector and matrix 
quantities needed to interpert the system in three 
dimensional space, with successive Figs. 6a - e 
showing different components of the measurement 
ma t r i x ; 

Fig. 7 is an equivalent schematic representation model of 
the application system of Fig. 1, including those 
model elements relevant to operation of the preferred 
embodiment; and 

Fig. 8 is a diagram similar to Fig. 7, but with additional 
details allowing for estimation of the accuracy of 
correction . 
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supporting table 8 resting upon the floor. Calculations 
performed within 5, transform force and torque measure- 
ments into touch location coordinates, which may be passed 
to a- controlling computer via cable 9, The controlling 
computer may then respond, via cable 10, with some altera- 
tion of the display image, such as brightening the image 
of a pictured push button whose coordinates have been 
touched . 

Before proceeding to a discussion of the features of 
the invention, a summary orientation of the force-measure- 
ment technique disclosed in my said copending application 
(also incorporated herein by reference) may be useful; 
through the improvement of the present invention is 
separate therefrom and useful also with other types of 
force-measurement systems as before explained. 

Returning, now, to the inertial interference 
suppression, cancellation or elimination concepts of the 
present invention, the stabilizer bar and bracket of Fig. 
1 was provided primarily to prevent overt motion within 
the tilt-swivel base 3, such that once the user has set 
the desired orientation of the monitor, there will be a 
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fixed spatial relationship between the face of CRT 11 and 
the platform device 5, 

The bulk of the supported mass resides in the rigid 
CRT 11 and its attached yolk assembly 12. Additional 
smaller mass concentrations are found in the monitor elec- 
tronics and power supply 13, and platform top plate 6, and 
a small amount is distributed throughout the plastic hous- 
ing and base. Significant elastic compliance is found in 
monitor housing 2, the platform device 5 itself, and 
especially in the tilt swivel base 3. Stabilizer 4 pro- 
vides some useful stiffness in parallel with 3, but not 
enough to reduce motions to insignificance. Given a 
perfectly stable bottom plate 7, the predominant two 
eigenmodes comprise rotations of the supported mass about 
horizontal axes passing slightly below the center of 3; 
one typically passing fore-aft, and the other left-right. 
The associated frequencies are close but not identical, 
and fall near 3 Hz for large, softly mounted monitors. As 
the relative importance of mass varies with the square of 
the distance f rom the axis $ the CRT and yolk assembly 
strongly dominate these modes. Numerous other eigenmodes 
reside at higher frequencies, some associated with motions 
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accurately modeled by a single, dominant, rigid mass, and 
some not; but all high enough in the application of this 
embodiment for linear filtration to suffice for removal of 
their disturbing effects. 

With acceleration readings available and properly 
exploited for motions of bottom plate 7 (or top plate 6), 
it has been found that the corrections relevant to the 
internal motions of the monitor/platform device system are 
mathematically identical to those for a motionless bottom 
plate as later explained. In the abscence of accelero- 
meters, however, motions of the bottom plate 7 and table 8 
must be explicitly considered. Folding tables typically 
allow large side-to-side shearing motion. Some varieties 
allow equally great motion front- to-back . Shelves and 
computer tops often allow large front- to-back nodding 
motion, due to torsion of the support surface. Resonant 
frequencies run 2Hz to 4Hz and are very easily excited, as 
are non-resonant motions during ongoing disturbance. 

Fig 2a depicts the major elements of a preferred 
accelerometer of suitable design for the particular pur- 
poses of the present invention, which is shown in exploded 



-19- 



assembly in Fig. 2b. A beri Ilium-copper foil 20 is etched 
or punched to the retangular outline shown, with mounting 
holes 21 and slits 22 formed in the same operation. 
Spring action resulting from slits 22 provides region 23 
with mobility perpendicular to the foil plane. Tabs 24, 
thermally isolated by slits 25, provide access for 
soldered or welded electrical connection. A single-sided 
PC board 26 comprises an insulating laminate 27 bearing 
conductive foil region 28, which serves as a capacitive 
counter electrode to region 23 of the free-standing foil 
20. Contact areas 29 are accessible along a thin strip 
not covered by foil 20 after assembly. Elements 20 and 26 
are held in opposition by spacing washers 30 and tubular 
rivets 31, with overtravel protection provided by a 
punched laminate bar 32. Wires 33 provide connection to 
the measuring electronics, and are typically routed some- 
what apart, so as to minimize stray capacitance there be- 
tween. The measuring system chosen, however, is largely 
unaffected by stray capacitance to the grounded enclo- 
sure. The dimensions and thickness of foil 20 are chosen 
to provide a resonance of 23 Hz. This is high enough to 
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provide good fidelity in rendering acceleration below 5 
Hz, yet low enough to avoid excessive aliasing energy 
within 5 Hz of the sampling rate. 

The resonant frequency is associated with a particu- 
lar distance such that region 23 sags under its own weight 
when the assembled acceleromet er lies horizontally. The 
thickness of washers 30 is chosen to be about half again 
this distance. The area of opposed conductors 23 and 28 
is chosen to give an operating capacitance range directly 
compatible with that of the main sensors (55 opposing 56 
in the upper and lower force-sensing platform castings of 
Fig. 3, later more fully described). The application at 
hand requires acceleromet ers with response free of signi- 
ficant amplitude or phase error from about 0.1 Hz to about 
5 Hz. The DC response characteristic of the design chosen 
for this embodiment is very helpful in this regard. 
Although DC coupling and capacitive technique are con- 
venient here, it should be noted that the method of the 

invention may employ acceleration responsive sensors of a 
wide variety of other types and principles of operation, 

if desired. 
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Fig. 3 depicts the force-sensing platform 5 seperated 
into its top plate casting 40, and bottom plate casting 
41, each with its associated sub-assemblies. Eight 
springs 42 have pressfit upper and lower still attachments 
43 and 44 respectively. Upper attachments 43 are held to 
bosses 45 by self- threading screws (not shown), with 
screws 46 (one shown) passing through lower casting bosses 
47 into threaded lower attachments 44. When thus 
assembled, and with connector 51 attached to circuit board 
53 through mating connector 52, the platform top plate and 
monitor are supported entirely through springs 42, the 
stiffness of wires 54 and additional supported monitor 
cables, such being insignificant compared to that of 
springs 42. Lower flange 48 surrounds upper flange 49, 
with lip 50 covering what would otherwise be an upward 
opening gap. This arrangement prevents incidental contact 
with keyboard, cables, or other small tabletop objects 
from passing extraneous forces to the top plate 40. Capa- 
citive sensors of displacement (and thereby of force) are 
located in pairs in the corners, each comprising the 
before-identified lower plate 55 which opposes an upper 
plate 56 upon assembly. The capacitor sensor plates 55 
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aud 56 comprise conductive foils on insulating subtrates. 
These sixteen substrates, plus the six acceler ometer' 
modules 57, are preferably adhesively mounted to 45-degree 
casting bosses, not shown. 



elements of Fig, 3 that pertain to the operation of the 
present invention. These are mounted in pairs, as shown, 
with axes of sensitivity inclined 45 degrees to the verti- 
cal, and each 90 degrees to the opposite member of its 
pair. The pairs are positioned to approximate the ver- 
tices of an equilateral triangle, and are turned to make 
the axes of sensitivity lie perpendicular to a line drawn 
from the triangle center. In practice, minor adjustments 
to this pattern may be made to avoid interferences. The 
triangle is made as large as may be accomodated within the 
base, to acheive a balance relationship between the rota- 
tional and translational sensitivities of the ensemble. 
This arrangement of accelerome t er positions and orienta- 
tions is, of course, only one of many, which will be evi- 
dent to one skilled in the art, that share the advantage 
of encoding the six degrees of freedom of rigid accelera- 
tion with excellent linear independence. 



The six accelerome ter modules 




are the important 
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A simplified schematic of suitable electronic conver- 
sion and calculating means for effectuating the functions 
of the invention is presented in Fig. 4, the circuits 
being carried on circuit board 53 within the platform 
device 5. This is similar to that used in my copending 
application and need not therefore be described in great 
detail. Circuits 85 and the multiplexer inputs 94 have, 
however, been increased in number (to fourteen), to acco- 
modate the six accele r ome t er channels 94b and eight dis- 
placement/force channels 94a. Timing circuit 83 has been 
modified to provide somewhat higher frequencies, suchthat 
forty complete conversions of all fourteen channels may be 
completed every second. Finally, microprocessor system 80 
has been enhanced with firmware to develop t irae d e r i va t i ve 
correction channels, and to employ acceleromet er and time 
derivative correction channels to improve force and touch 
location data made available over communications cable 9, 
Fig. 1. The capacitance 57 in circuit 85, furthermore, 
will be replaced by 55 and 56 for the inputs 94a. 

It is now in order to explain how circuits and compo- 
nents operate to obtain the desired functions. In Fig. 5, 
a flowchart of firmware is presented that may be used in 
the microprocessor system 80 to transform the sensor input 
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data into force and the torque substantially free of iner- 
tial error. The subsequent processing to yield touch 
location coordinates, as described in my said copending 
application, is not germaine to this operation and there- 
fore, .not repeated in this application. 

The conversion methods of Fig. 4 sample the fourteen 
input channels at fourteen different points in time, 
more or-less evenly spaced throughout the interval of one 
full sample cycle. Within each fourteenth of a cycle, the 
exact moment of measurement depends somewhat on the actual 
value measured. The resulting differential group delay is 
removed by processing step 60. Step 60 may comprise four- 
teen different two-tap finite impulse response ("F.I.R.") 
filters, each performing a linear interpolation of the two 
most recent component values to a common equivalent sample 
time. At the cost of slightly greater total delay, 
three-tap quadratic interpolating filters may be used to 
acheive substantially less differential amplitude 
response • 

In step 61, Fig. 5, by contrast with step 60, the 
same response is employed in each of the fourteen filters. 
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Th i s response is designed to remove frequencies above 
about 5 Hz. In one embodiment, for example, 12 tap 
coefficients form a triangular impulse response, compris- 
ing six values in descending sequence. This is equivalent 
to successive applications of a six-wide and a seven-wide 
square impulse response filter, which give notches at, for 
example, 5.7 Hz, 6.7 Hz, and integral multiples of these. 
As noted in my said copending application, a linear filter 
which is applied uniformly to all components of the 
applied force and torque vector leaves the computed loca- 
tion of force application unaffected, except for possible 
delay in response. Further, the linear nature of the pro- 
cessing steps in Fig. 5 provide that the uniform applica- 
tion of a given filtering characteristic to all channels 
in step 61 is equivalent to applying the same characteris- 
tic instead to the six channels output from 70, insofar as 
final output 74 is concerned. 

The compensatory differential filtering of step 60 is 
thus provided first to assure that the effective aggregate 
filtering of all channels is identical, and then an add- 
itional common characteristic is introduce in step 61 to 
remove frequencies higher than needed for adequate re- 
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sponse speed, which frequencies may also be too high for 
accurate handling under certain simplifying assumptions 
exploited by later processing steps (such as the use of 
second central difference to approximate the second deri- 
vative . ) 

Following step 61, the eight channels responsive to 
the force sensors are standardized by linear transforma- 
tion 62, and the six channels responsive to the accelero- 
meters are standardized by linear transformation 63. 
These transformations are assigned the task of compensat- 
ing for mechanical and electronic component variation from 
one particular platform device to the next, so that if 
multiple devices are used in a set of identically con- 
figured installations, the installation specific calibra- 
tion data may be identical across the set. Transformation 
62 has an additional task of reducing the eight force sen- 
sor channels to six, in such a manner as to eliminate the 
effect on any of these channels of distortion of the 
mechanical structure of the platform device. The forty 
eight coefficients of combination required by 62, and the 
thirty six required by 63, are determined by a factory 
calibration procedure and stored in a non-volatile memory 
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which is part of the microprocessor system 80. The matrix 
of coefficients for 62 may be described as in my said 
copending application. Determination of the matrix of 
coefficients for 63 is simplified by the fact that the 
only significant source of variation resides in the over- 
all scalar sensitivity of each channel, as opposed to the 
position and direction of its axis of sensitivity. These 
channel sensitivities may then be determined by supporting 
the platform device briefly at each of two fixed angles, 
and noting the differing effect of the acceleration of 
gravity on each channel for the two cases. A diagonal 
matrix comprising the ratios of the expected differences 
to the observed differences may then be used for 63. 
Alternatively, 63 may be the product of this diagonal 
matrix with a fixed matrix describing a transformation to 
familar rectangular coordinates. 

The six components of the most recent previous output 
of 62 are held in delay 64, while the second most recent 
are held in delay 65. Combiner 66 sums, component by com- 
ponent, the current and doubly delayed outputs of 62, 
while subtracting twice the singly delayed output, yield- 
ing second central difference signal 72. Delay 67 is pro- 
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vided to match the group delay of forming the second cen- 
tral difference, such that the signals 72 and 73 remain 
accurately aligned in time. The second central 
difference, before mentioned, is employed as a good low- 
to-medium frequency approximation to the true second deri- 
vative, with the advantage of small group delay. Five or 
seven point approximations may be used in other embodi- 
ments for accuracy to higher frequency, but at the cost of 
greater response delay. The uncorrected force readings, 
available with the appropriate degree of delay from 64, 
then pass to combiner 70, wherein corrections from trans- 
formations 68 and 69 remove the inertial interference sig- 
nal. The coefficients of combination 68 and 69 are deter- 
mined in an installation dependent manner by a user cali- 
bration procedure, performed in conjunction with a soft- 
ware program executed on the controlling computer with 
which the platform device communicates. At the end of 
this procedure, coefficients calculated within the con- 
trolling computer are passed to microproccesor system 80 
for storage in its non-volatile memory. 
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The touch locating application for which this embodi- 
ment is intended involves a series of short term measure- 
it 

ments of changing forces, separated by periods of quie- 
scence. An "auto-zeroing" feature 75 can therefore pro- 
vide economy by relaxing the requirement for long term 
stability against measurement offset drift. A variety of 
heuristics may be employed for determining quiescence, 
such as the failure of a number of most recent samples to 
contain any components with a range of variation exceeding 
a certain minumum. Quiescent samples identified in this 
way (all of whose components shouJLd nominally be zero) are 
accumulated into a set of baseline values. These baseline 
values are presented continuously to combiner 70 for sub- 
traction of static offset from outputs 74. By accumulat- 
ing only a fraction of each quiescent sample into the 
baseline values, a simple recursive filtering effect may 
be obtained which prevents the auto-zeroing logic 75 from 
degrading the signal-to-noise ratio of outputs 74. 
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CALIBRATION PROCEDURE 

The calibration procedure for transformations 68 
and 69, Fig. 5, is performed as follows: 

- At installation time, the user loads and executes 
a software program supplied with platform device 5 on the 
controlling computer C to which it is attached. This 
program, so operating, is able to pass operating 
instructions and action acknowledgements to the user 
through display device 2. The initial steps of this 
program include establishing communication with platform 
device 5 through cable 9, Fig. 1, and instructing the 
device to enter a calibration mode in which all 
components of signals 71, 72, and 73, Fig. 5, are 
streamed to the controlling computer, rather than the 
usual touch location information. The program then 
monitors the three translational components of signal 71 
for a short while, to establish their quiescent zero 
points. The squared magnitude of applied translational 
force is thereafter obtained from the sum of the squares 
of the deviations of these three components from their 
zero point levels. 

The program next leads the user through a series 
of four to six mechanical disturbances of the system, 
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ciuring which it gathers a matrix of twenty four to thirty 
six measurement vectors, of eighteen components each, 
from which transformations 68 and 69 will be calculated. 
For the first four disturbances, the user 1 is instructed 
to apply finger pressure P at certain approximate 
locations on display monitor 2, Fig. 1. In each case, 
upon detecting an applied translational force of 
appropriate magnitude and stability, the user is prompted 
to remove his finger abruptly, and keep it off until 
instructed to proceed to the next location. During the 
first second after release, six measurement vectors are 
taken from the forty or so available and are retained for 
later processing. The interval between successive 
retained measurements is increased from five to nine 
samples over the interval, to eliminate the chance, for 
some resonance, of strobing only zero crossings. 
Simultaneously, a running sum is kept for the sample 
values of each component during the interval. After the 
six measurement vectors have been gathered, they are each 
diminished by the average vector obtained by dividing the 
running sum vector by the number of samples summed. 

After further fluctuations in the uncorrected 
translational force samples have died away to a low 
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level, the user is instructed to continue with the next 
location. The four locations and directions of finger 
pressure prior to release are typically as follows: the 
top right corner of the monitor 2 pressing toward the 
rear; the top left corner pressing toward the rear; the 
top right corner, around to the side pressing to the left; 
and the bottom right corner, around to the side pressing 
to the left. However , neither the number nor the 
positions of the locations actually used is critical. 
What is needed is a sufficient number of measurements of 
non-zero free motion of the system, in which each easily 
excited eigenmode (or more precisely, conjugate eigenmode 
pair) makes a number of contributions to different 
measurement vectors at randomly differing amplitudes and 
phases. In addition, the number of data elements 
gathered from correction channels (components of 72 and 
73, Fig. 5) should substantially exceed the number of 
coefficients to be determined, so that full constraint of 
the latter may be insured, despite the stochastic nature 
of the independence of eigenmode representation. 

In addition to the twenty four free-motion vectors 
so gathered, an attempt is made to add two more groups 
corresponding to motions driven from the support side. 
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The user is instructed to shake the table gently from 
left to right, A group of six measurement vectors is 
collected, as above, except that collection commences as 
soon as the program detects a sufficient level of 
disturbance, and the user is instructed to continue 
shaking through the measurement interval* If sufficient 
disturbance is unavailable with gentle shaking, the 
program will offer no immediate acknowledgement — a 
condition to which the user has been instructed to 
respond by pressing a key on the controlling computer's 
keyboard. Finally, instructions are given to repeat this 
process by shaking from front to back. The driven 
measurements collected in this manner show response to 
frequencies not constrained to the eigenf requencies of 
free motion; they may be particularly helpful in 
constraining the coefficients of acceleration correction 
matrix 69. 

Were the support soft enough to allow slow, 
substantial motion in all six degrees, the associated 
measurement vectors would together tightly constrain all 
aspects of transformation 69. Few tables of any kind, 
however, allow significant motion in vertical 
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translation ♦ Similarly, the "visible" (i.e., it affects 
the measurements) excitation of all eigenmodes 
corresponding to the principal supported mass would fully 
constrain transformation 68. Yet the presence of filter 
61 virtually assures that at least some of the higher 
frequency modes will be invisible in the data. 
Nonetheless, effective coefficients for transformations 68 
and 69 can be found from the measurements vectors 
collected as described above, in the sense that they 
correct well for the disturbances actually occurring in 
the installation at hand, while manifestly having no need 
to remove disturbances never found in the data. 

After completion of the measurement-collecting 
phase, the calibration program enters a calculating 
phase, during which it finds such coefficients for 68 and 
69. These are then transmitted to platform device 5 for 
storage in the nonvolatile memory of microprocessor 
system 80. 

At this point, the program may proceed to 
determine static calibration coefficients necessary for 
touch force location, as described in my said copending 
application. 
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Consider the equation 

f=f-Af-Ba, (1 
c u u 

*e "f " and "f " denote force measurement six-vectors 
c u 

ind 71, corrected and uncorrected for inertial errors, 

>ectiiirely; where "a" denotes acceleration six-vector 

whers* "A" and "B" denote the six by six constant 

:icesi;~68 and 69, and where the number of dots placed 

c an Instance of a variable quantity denote the order 

lime- -derivative with which the variable is taken in 

instilnce. Equation 1, then, closely approximates the 

culat$;ons of Fig- 5 subject to the bandwidth 

itatign intentionally placed on all variables. 

It has been determined empirically that matrices 

and "B" may be found such that "f " very much more 

c 

sely approximates the true unknown applied force "f" 

n does "f • It is also demonstrated in the 
u 

einafter presented analysis, that this is a reasonable 

>ectation. Where, as for the calibration measurement 

:tors, the applied force "f" and the desired estimate 

it "f " may both be taken to be zero, we have the 
c 
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but avoids possible problems of overflow due to 
unexpectedly large coefficients- Additionally, should 
modes of disturbance poorly excited during calibration 
ever be strongly excited in practice, the less 
sophisticated determination of "C" may greatly enlarge 
the error, where the alternative at least leaves it 
alone . 




THEORY UNDERLYING THE INVENTION 

To understand how the method of the invention may 
be varied in embodiment, and how it may be applied to 
differing force measurement problems, it is necessary to 
examine the underlying theory as part of the teaching. To 
do this, we will examine how lumped parameter mechanical 
models can be used to analyze the touch location 
application, both with regards to the correct functional 
form for inertial correction equations, and with regard to 
the limits on the accuracy improvement which may be 
expected. 

Mechanical schematic diagrams are usefully 
employed, such as Fig. 7^/ Fig. 8^~ and - F i g. —9a. While 
correct equation forms may be readily extracted from 
these by viewing them as one-dimensional physical 
diagrams, they become full six degree representations of 
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the three-dimensional problem with correct interpretation 
of the symbols and associated quantities* Capital 
letters denote six by six real, symmetric, constant 
matrices, while lower case letters denote variable 
six-vectors containing both translat ional and rotational 
information. The matrices represent not only all 
degrees of those parameters normally associated with the 
component elements (e.g. spring stiffness, mass, moments 
of inertia, magnitude of damping) , but also aspects of 
relative component position, such that in conjunction with 
these matrices , the diagrams need convey only the 
topology of interconnection to fully specify the 
problem. 

Fig. 6 assists interpretation of vector-valued 
variables and model parameters in three-dimensional 
space. The choice of reference system for these 
quantities is arbitrary, in the sense that the forms of 
the mechanical circuit node equations, and results derived 
from them, are independent of the reference system used. 
Any analysis which requires numerical estimates of the 
component values, however, must choose a convenient 
reference system, and then apply that same one 
consistently to all quantities in the problem. The 
required reference comprises a set of coordinate axes 
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100, Fig. 6a, together with a fixed point R in space. 
Note that axes 100 define not only orthogonal 
translational degrees "1", "2", and "3", but rotational 
degrees "4", "5", and "6", the latter three being the 
outward looking clockwise rotations about the first 
three, in order. This enumeration order is depicted at 
100 as applied to the components of force "f", but 
applies analogously to all six-vector quantities. 

In Fig. 6a, the set of forces 101 acting on rigid 
body 102 is summarized as to effect, in total state of 
force six-vector "f". An agent at R, supported by rigid 
extension 103 from solid support 104, reports in terms of 
coordinate axes 100, the force and torque components 
required to maintain rigid extension 105 of body 102 in 
equilibrium. The negative of these components, 
enumerated in the order shown at 100, constitute the 
vector "f". (Reversing sign reports the force felt, 
rather than the one applied.) Note that the value of 
"f", while strongly dependent upon the position chosen 
for R, is independent of the path of 105 or its point of 
attachment to 102. Note also that many different force 
patterns 101 may lead to the same value of "f", but that 
the effect of any such force pattern on the dynamics of 
body 102 is completely determined by the value of "f" . 



Fig, 6b shows rigid body 102 displaced from its 

equilibrium position "x " by displacement six-vector n x", 

e 

All possible rigid extensions of body 102 define a mobile 

but rigid auxilliary frame of reference termed a frame of 

rigid subsystem 102; or, more briefly, the frame of 102. 

At R, the frame of 102 has been marked with a copy 10 6 of 

coordinate system 100. After displacement, the first 

three components of "x" may be read as the displacement 

of the marked copy of "R" with respect to "R" itself. 

Some care is required in addressing the representation of 

rotations; but in the limit of small displacements, the 

final three components may be read as that rotation of 

the original coordinates at "R", about each axis in turn, 

which brings the other two axes closest to their marked 

copies. Thus the best aligning rotation about the third 

axis develops the component "x ", as shown. 

6 

Alternatively, a rotation vector (not shown) may be 
defined which points parallel to the axis of rotation in 
the clockwise direction, with length proportional to 
angle of rotation about this axis. Components "x 

4 

"x and "x " are then the projections (not shown) of 
5 6 

this vector on the first, second, and third coordinate 



axes, respectively. Importantly, in the limit of small 
displacements, the sum (component by component) of two' 
displacement six-vectors accurately reflects the 
compounded displacements. 

In Figs* 6c, 6d, and 6e, a mechanical system is 
shown composed of rigid and potentially mass bearing 
subsystems 102, 107, and 108, coupled by massless elastic 
element 109, with characterizing matrix S, and massless 
viscous dissipation element 110, with characterizing 
matrix D. For brevity, elements such as 109 are termed 
springs, and elements such as 110 are termed dashpots; 
these terms imply nothing, however, as to the physical 
form of the elements in question. The spring symbol at 
109, both here and in the schematic diagrams to follow, 
is taken to denote the presence of stiffness in all six 
degrees, not just along the axis of the symbol. Also, it 
may be a summary representation of the effects of a 
number of physically distinct elastic bodies which 
interconnect the rigid subsystems (e.g. 102 and 108) to 
which it attaches. Analogous comments apply to the 
dashpot symbol at 110. Note that while later diagrams 
are purely schematic, Fig. 6 is partly geometric for 
expository reasons. In Figures 6c, 6d, and 6e, spring 
109, dashpot 110, and mass 102 are dissected free of the 



rest of the system, and connected with massless rigid 
extensions for measurement of their matrix parameters 
"S", "D" , and "M" as seen at R. 

Returning to Fig. 6c, the agent at R imposes small 
static translations or rotations of 111 with respect to 
112, in each case along one of the six degrees of 
coordinates 100, while at the same time applying whatever 
forces and torques are necessary to preclude motion along 
any of the other five. In each case, the force six-vector 
felt (or change in force six-vector felt, if the initial 
position is not relaxed) , divided by the magnitude of the 
displacement, is entered as that column of the six by six 
matrix "S" corresponding to the number of the component 
displaced. It follows from this definition, then, that 
if "x" is a displacement six-vector of 111 with respect to 
112, from an initial equilibrium position, then the force 
six-vector felt through 111 is given by the product 
"S x". Since "S" is necessarily symmetric, this is also 
the product "x S". A force "-S x" is felt through 112. 

In Fig. 6d, similarly, the agent at R imposes 
three translational velocities and three rates of 
rotation upon 113 with respect to 114. In Fig. 6e, the 
agent at R imposes three translational accelerations and 
three rotational accelerations upon 115 with respect to 



116. In both cases, the motions are constrained to one 
coordinate degree at a time, as in Fig. 6c. As before, 
coefficients for the matrix columns of "D" may be 
determined in order, by taking the force six-vector felt 
for each motion (the negative of that applied) , in order 
of degree, and dividing its components by the magnitude 
of the velocity applied. Coefficients for the matrix 
columns "M" may be determined in order, by taking the 
applied force six-vector for each motion, in order of 
degree, and dividing its components by the magnitude of 
the acceleration applied. Like "S"> "D" and "M" are also 
always symmetric matrices. 

"S" and "D" are defined through "felt" forces, 
while "M" is defined through the forces "applied". This 
leads to a sign convention in which the diagonal 
components of "S" and "D" are negative, but also one in 
which terms in node equations of the general form 
"f = M a' 1 are all conceptually of positive sign. Actual 
negative signs appear only in writing terms for which 
the available displacement variable is used in a 
negative sense; that is, it describes the motion of some 
other node with respect to the node of interest. 

In the intact and excited mechanical system, let 
"x" be the six-vector displacement from static 



equilibrium of 108 with respect to 102, let "y" be the 
six-vector displacement of 102 with respect to 107, ancl 
let "z" be the six-vector displacement of 102 with 
respect to the inertial frame of R, 104. Then all of the 
forces acting on rigid subsystem 102 are the force 
six-vector exerted by 109 on 102, which has the value 
"-S x", and the force six-vector exerted by 110 

on 102 f which has the value "D y" . The total force 
acting on 102 necessary to explain the observed 

acceleration of 102 is "M z " . 

Thus the node equation for rigid subsystem 102 
may be written: 

• * • 

-Sx + Dy = Mz (7 
In similar manner, equations for all the nodes of any 
lumped parameter system may be written. These may be 
solved to provide a description of the system dynamics, 
wherein combinations of force six-vectors of interest and 
their time derivatives of various orders are equated to 
combinations of dispacement six-vectors of interest and 
their time derivatives of various orders. 

This completes the analogy between the six-vector 
desciption of six-degree per element mechanical systems 
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and the usual scalar valued interpretation of the 
equations of linear system theory, allowing the extensive 
results of the latter to be brought to bear on the more 
compact six-vector notation. {In scalar form, the 
equations are approximately thirty six times as 
voluminous . ) 

It is now appropriate to turn to Fig. 7, which 

depicts a lumped parameter model which approximates the 

dynamics of the system of Fig. 1. Platform base 7, the 

top of table 8, and any other objects supported by 8 are 

considered to form a rigid subsystem which is node 134 

(with mass matrix n M ") . Display monitor 2, exclusive of 

1 

its tilt-swivel base, is considered to form a rigid 

subsystem which is node 131 (with matrix "M ") . The legs 

2 

of table 8 are considered to be an elastic coupling 135 

(with spring matrix "S ") between node 134 and solid 

1 

floor 136. Elastic coupling 133 (with matrix "S ") 

2 

represents the joint effect of the eight springs 42. 

Elastic coupling 132 (with matrix "S ") represents the 

3 

joint effect of tilt-swivel base 3 and stabilizer bar 4 



of Fig. 1* Platform top 6 is considered to form a rigid 

subsystem which is node 139. This is taken for the time 

being as massless, as are bodies forming elastic couplings 

132, 1'33 and 135. The force "f" to be measured acts upon 

monitor node 131, while interfering force "f ,f acts on 

t 

tabletop subsystem 134. The available data comprise the 
displacement six-vector "x" of subsystem 139 with respect 
to subsystem 134, and optionally, acceleration six-vector 
"a" of subsystem 134. In most actual embodiments, of 
course, available data are invertible memoryless 
transformations of these. 

Massless node 139 may be described by equation 
140, Fig. 7, which leads immediately to expression 141 
for the six-vector displacement of node 131 with respect 



to node 134. ^Equations F±^ 7b and 7c of nodes 134 and V" 
131, respectively, may then be written without reference 
to displacement "z\ Also, the acceleration "a" measured 
for node 134 must correspond to model values as shown in 
Trtb^ 7d. 

Manipulation of Fig^, 7b and 7c yields F±q^ 7e. 
, This equates a function of unknown forces only, with a ^ 
function of known displacements only. Combining Fiq. 7c V 
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and 7d, however, gives ^i^. 7f. This latter equation has 

simplified dependencies which provide some of the 

advantages of incorporating accelerometers in accordance 

with the techniques of the invention. 

Since the uncorrected force readings "f " 

u 

correspond directly to "-S x", ^"ig. 7f may be recast in 

2 

terms of "f " as seen in ^^ig^ 7g, ^Mcf. 7g corresponds 
u 

directly to the form of previously presented equation 1, 
thus explaining the excellent results obtained in 
practice in its application to the system of Fig, 1. 

In these equations, the ratio of a second 
derivative term to the zeroth order term of the same 
variable (or of fourth order to second) is always seen to 
give a "coefficient ratio matrix" with the general form 

-1 -1 

M M S The quantity "M S " describes the free motions 

of a system consisting of nothing but the mass and spring 
referred to, the far end of the spring being clamped . The 
eigenvalues of this matrix are the squared radian periods 
of the resonances of the mass so supported. Where the 
matrix is characterized by a single frequency "v", the 

-1 . . 

terms "x" and "M S x have comparable magnitude when 



"x" is dominated by components with frequencies near 
"v" . For a frequency of one-tenth "v", however, the 
zeroth order term dominates the second order term by a 
factor- of one hundred (the second derivative scaling as 
the square of frequency) . In general, though, the matrix 
has a broad range of eigenfrequencies, and many 
appications, such as touch location, have significant 
sensitivity to error from excitation of all modes. It 
may be concluded as a general rule, then, that the lower 
order term can be expected to dominate the term of order 
two higher by at least the square of the ratio of the 
lowest eigenfrequency of the coefficient ratio matrix to 
the frequency of interest . As an indication of the 
relative error in one degree of freedom of the final 
system output occasioned by omission of a higher order 
term, this ratio may prove pessimistic by a modest 
factor, since the system is usually designed for good 
sensitivity in the sense that typical "x" values have 
large projections onto that degree' s axis in the vector 
space of "x", while the error components project 
more-or-less randomly. 

The ratio of uncorrected force terms in >ig. 7g 
(or of displacement terms in £<g. 7f) describes the 
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resonances of the mass of 131, Fig. 7 , with the series 

connection of springs 132 and 133. Interestingly, the'se 

are just the free motions of display 2 in Fig, 1 when the 

table , is absolutely solid. The accelerometers, as used in 

the present invention, thus have created a virtual solid 

platform; "f "S ", and "M " have entirely disappeared 
t 1 1 

from view. The acceleration term makes the system immune 

to inertial interference from the support side, while the 

second order term makes it immune to interference from 

motion induced by changes in the unknown force. 

In a configuration where the lowest of the 
resonances just described is sufficiently high, the second 
order term may be ommited, and the acceleration term 
retained to deal with table tilt and motion. When the 
table in reality is solid, of course, the acceleration 
term (and the accelerometers) may be omitted, and the 
second order term retained to deal with a heavy monitor 
on a soft base. 

The absence from iN^g. 7g of constants pertaining 
to the mechanics of the table means that calibration 
matrices exist to correct a given installation for all 
possible support surfaces, the particularities of the 
installation being limited to the mechanics of the 



display and its relationship to the platform device. To 
achieve such a "portable" calibration, however, the ■ 
procedure described earlier must be performed on a 
surface mobile enough to allow acceleration in all 
degrees of freedom that might ever plausibly be seen. 



without accelerometers, full correction can be obtained 
from information impicit in the history of the 
uncorrected displacement (equivalently , force) data 
itself, subject to certain restrictions. The most 
important limitation resides in the absence of any ability 
to remove error resulting from motion forced through the 
support. Thus when the table is not very solid, the 
frequencies associated with coefficient 143 are quite 
low, and table forces "f " must be kept very small. 



Also, these equations require 135 to be a true linear 
elastic member, which can be a questionable assumption, 
due to stiction with, or bumping against, adjacent walls 
and furniture (the accelerometers mask this assumption) . 
Again, when the table is not very solid, the frequencies 
associated with coefficient 142 are low; thus corrections 
from the fourth order derivatives must be made and must 
be calibrated (which may be done in the same manner as 



Returning to 




7e, it can be seen that even 



t 



previously described for the other correction channels) . 
Thus also, the calculating means of microprocessor system 
80 must be modified to integrate out the second order 
force term on the left. This requires calibration of 
coefficient 142. This calibration is still preferably 
performed by in situ, empirical means, but may require 
additional steps with known non-zero forces, which may 
require equipment. 

Further refinements and accuracy in achieving the 

immunity to inertial interference effects in accordance 

with the invention may be attained with the aid of the 

somewhat more realistic lumped parameter model of the 

dynamics of the actual system of Fig. 1 represented in 

Fig. 8. Node 131 of Fig. 7 has been split into 

independently moving masses 150 and 151. Node 151 (with 

mass matrix "M ") represents the CRT 11 and yoke assembly 
4 

12 (and such parts of the frame and housing as are best 

treated as co-moving with 151) . Node 150 (with mass 

matrix "M ") represents the power supply and electronics 
3 

13 (along with such parts of frame, housing, stand 3 and 

stabilizer 4 as are best treated as co-moving with 150) . 

Spring 153 (with spring matrix "S ") has been introduced 

5 




to describe elastic deformation of the frame and housing 

of the CRT display between 11 and 13. Spring 132 

(Fig. 7) has been replaced by spring 152 (with spring 

matrix- "S ") in recognition of the fact that spring 132 
4 

of Fig. 7 implicitly incorporated some of the effects of 
spring 153. Accelerometer measurements are presumed to 
be available. , . 

The node equations , presented if*-£-ig% 8b and 8c, 
lead to force equation F^j. 8d: ~ ^~*wy>- ffy *\ { "\ s 

Two conclusions are drawn from Fa>g. 8d. One is 
that in some applications and embodiments, corrections of 
greater accuracy may be achieved within the method of the 
invention by including corrections from higher order 
derivatives of the uncorrected force measurements, and/or 
from derivatives of the accelerometer measurements. In 
this instance, calculation and calibration may be 
performed exactly as previously described, except for the 
necessary widening of matrix "C" and heightening of "R 

cc 

The second conclusion is that the errors of the 
Fig. 7 model corrected in Fig. 8 are quite minor, and 
that for the actual application of Fig. 1, a scheme more 
elaborate than that of Fig^ 7g is probably not 



justified. First, consider the sub-terms 159 and 160 
(Fy£g. 8d) that comprise the coefficient of second order 
term 157. The ratio of coefficient sub-term 159 to the 
coefficient of term 158 describes the resonances of node 
151 with the series connection of springs 133, 152, and 
153. The same ratio for coefficient sub-term 160 
describes the resonances of node 150 with 133 and 152 in 
series. The lowest frequencies associated with 159, 
then, are much lower than those of 160; partly because 
the scalar mass of 150 is typically at least several 
times smaller than that of 151, and partly because the 
center of gravity of 150 is about three times closer to 
the center of rotation of its effective spring. Lowest 
frequency modes in both cases are primarily rotations 
about the spring centers, and depend upon moments of 
inertia that scale with the square of the center of 
gravity to spring center distance. The coefficient of 
157 is thus dominated by 159, and the ratio of the 
coefficient of 156 to this is then close to expression 
161. Since the mass of 151 is substantially larger than 
that of 150, the ratio of the coefficient of 154 to that 
of 155 is also approximated by expression 161. These two 
ratios determine the importance of terms 154 and 156, the 
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two terms that are outside the original form of equation 
1. The mass of the power supply, electronics, etc. is* 
relatively low, and the stiffness of the housing and 
frame -is relatively high, especially in the sense of 
there being no axes of spring rotation passing far from 
the center of gravity of 150 to create low frequency 
torsional modes. The lowest frequencies associated with 
expression 161 are therefore found to be in the range of 
tens of Hz. The principal frequencies involved in 
locating even very rapid touches are in the vicinity of 
two to three Hz for the preferred embodiment. Thus 
correction terms 154 and 156 are roughly one hundred fold 
smaller than those already allowed for in equation 1. 

The approach of Fig. 8 may also be used to 
evaluate the effect of bodies other than the power 
supply, to see if the failure of their mass to move 
rigidly with that of the CRT might require corrections 
outside the original form of equation 1. In the 
application of Fig. 1, the top plate 6 is the only other 
candidate for a mass of significant effect. Thus, Fig. 8 
may be reinterpreted with node 150 taken as the mass of 
top plate 6 (Fig. 1) instead of power supply 13. 153 is 
taken as base 3 with bar 4, 151 is taken as the bulk of 



display 2, and the equations are slightly simplified by 
-1 

setting "S " to zero. For the application of Fig, 1, 
4 

the mass of the top plate is about ten times less than 
that of the monitor. (The ratio may be somewhat less for 
a small monitor, but a small monitor suffers less 
interference and requires less accurate correction.) The 
top plate center of gravity is several times closer to 
the spring center of base 3 than is the center of gravity 
of the monitor. Together , these imply that the lowest 
resonance of expression 161, as here- interpreted, is 
about ten times higher than that of the monitor base 
combination. As this latter is known to be about three 
Hz, it can be seen that the correction terms for top 
plate mass would again affect three Hz. force information 
only about one one-hundredth as much as the correction 
terms present in the original form of equation 1. 

Dissipation features have been omitted from the 
model systems considered to this point, since typical Q's 
observed for the low frequency modes requiring the 
greatest correction in the application of the preferred 
embodiment are in the range of fifty to one hundred, 
viscous dissipation features, however, are readily 
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incorporated by following the general approach of Figs, 7 

and 8, This introduces terms with derivatives of odd ■ 

order into the force equation derived. These terms have 

coefficient ratios to the next lowest order terms that 

express first order relaxation time constants; the 

reciprocal of the longest in each case being the radian 

frequency at which the higher term may acquire comparable 

importance to the lower. 

Progressively more elaborate models aquire terms 

with progressively higher derivatives, and force 

derivatives are likely to appear on the left side. In 

some cases, the latter may be banished in a series 

solution formed by repeatedly differentiating and scaling 

the equation, and subtracting the result from itself; 

utility of such depending upon the frequency range of 

series convergence. F>g. 7f may indeed be extended by 

this means to incorporate a single dashpot parallel to 

132, as the interesting relaxation times (e.g. of the 

plastic tilt-swivel base taken by itself, when briefly 

deformed and released) correspond to frequencies well 

above the range to be corrected. The lowest frequencies 

of 142 (F^i^ 7e) , however, are much too low for the 

series solution of t*ig^ 7e ("f H assumed zero) to have a 

t 



useful range of convergence. 

In any case, extensions to the original form o'f 
equation 1 involve the addition of more terms for 
derivatives of some of the other orders of the variable 
six-vectors. These are generally of relatively low 
order, as term significance drops off rapidly with 
increasing order for the frequencies of interest. Series 
related terms need usually include only the base member, 
and the highest order terms arising in a model may be too 
small to need inclusion. 

When empirical calibration is used to determine 
the coefficients of equation 1, in original or extended 
form, best fit coefficients will be found given the terms 
available, regardless of what model or expectations 
motivated the inclusion of such terms. Thus Eig^. 7f and 
F^kcj. 8d provide somewhat different physical 
interpretations for the coefficient of the second 
derivative of the uncorrected force six-vector. 
Progressively more detailed models would provide 
progressively more elaborate and accurate 
interpretations. Careful spring and mass measurements, 
taken from progressively finer dissections of the system 



of Fig, 1, could in principle be used to calculate this 
progression of coefficients* Each such matrix in turn 1 
would have a somewhat different numerical value, but 
converging rather rapidly toward the single "true" value 
for the second order coefficient. Yet this latter is 
just the value found — given that the coefficients of 
other necessary terms are also being evaluated — by the 
empirical , in situ calibration method hereinbefore 
described. Determination of what other terms are 
necessary may be approached theoretically, as has been 
seen. It may also be approached experimentally, by the 
simple expedient of trying forms of equation 1 with 
differing sets of derivative orders present, to see what 
works best for the class of mechanical system relevant to 
a particular application. It should be noted that 
ommitting an unneeded term of lower order may be as 
important as including a necessary term of higher order, 
as ambiguities of solution can otherwise exist which 
preclude the convenience of a calibration in which the 
unknown forces applied to the excited system may all be 
zero » 

Finally, in the preferred embodiment of the 
invention, evaluation of the correction terms in 



equation 1 is subject to known inaccuracies: a slight 
high frequency roll-up in the accelerometer channels 1 
reflecting their twenty three Hz resonance; the 
variations in the filtering components of the circuits 91 
(Fig. 4); and the use of the second central difference 
for the second derivative. Each of these is a source of 
almost two percent error in the corrections. (Since 
these errors are independent their combined expected 
value is around three percent, consistent with the before 
mentioned factor of thirty improvement attained.) 
Embodiments seeking improvement of correction for this 
application thus begin with improving the accuracy with 
which the present correction terms are evaluated, rather 
than by adding new terms to equation 1. 

Therefore, in an alternative embodiment, sampling 
frequency may be increased, the accelerometer resonant 
frequency may be increased, and the two output capacitors 
of each circuit 91, preferably all in_ similar ratio, may 
be decreased. These changes improve the accuracy with 
which the correction terms of original equation 1 are 
evaluated, thereby further improving overall system 
rejection of inertial interference. 
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In still another embodiment, timing circuit signal 84 
in the circuit of Fig. 4 is forced to a somewhat different 
amplitude by factory test apparatus, and the step 
responses of circuits 94 monitored, such that correcting 
inverse responses may be convolved into the filters 60 of 
Fig* 5. A correcting inverse over the frequency range of 
interest (up to about ten Hz) for the accelerometer 
resonance is likewise convolved into the accelerometer 
channel filters of 60. The filters of 60 thus modified 
are convolved into the filters of 61. The second deriva- 
tive approximation may be expanded .to a five tap filter. 
The linearity of operations is exploited to change the 
order of differentiation and matrix multiplication, so 
that the second derivative filter may be convolved into a 
copy of the eight force responsive filters of 61, as modi- 
fied. Related matching delays may also be added to the 
other channels, and two copies of transformation 60 must 
be executed to provide both signal 71 and signal 72 of 
Fig. 5. Each of the resultant filters 61 is then 
shortened in group delay by truncating the (very small) 
lead coefficient. In this manner, correction accuracy is 
improved . 
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In either of these last embodiments, the low-pass 
cutoff frequency of filters 61 may be raised, allowing 
lower group delay. In this manner, improved correction 
may be exploited to provide faster system response. 

While the invention as applied to the touch-input 
computer and related CRT or similar displays employing 
touch force location measurement techniques has been illu- 
stratively described, preferably embodying each of the 
novel accelerometer , derivative correction and calibration 
techniques of the invention, there are instances where the 
combination of all such techniques may not be required for 
inertial interference suppression, as before mentioned. 
Concerning the independ ence of the accelerometer and 
derivative connection techniques, there are touch screen 
systems such as some monitor or desktop systems without 
control of the support surface where only the accelero- 
meter - based connection may be most useful. For a 
general purpose retrofit unit, on the other hand, a major 
problem resides in large or softly mounted monitors. Such 
units would benefit from the derivative - based correction 
only, at least for users with a reasonably solid desk or 
table. The novel calibration procedure is also more 
generally useful. 



From another approach or viewpoint, a feature of the 
invention resides in its novel form of an effective linear 
filtering, at least as applied to force measurement. 
While,, as earlier discussed, simple forms of so-called 
linear filtering have been proposed in other types of 
applications, such^unlike the present invention, generally 
attenuate all energy near the frequencies of distortion, 
rather than correcting the distortion. The present inven- 
tion, as previously explained, quite to the contrary, 
attenuates the error of the signal but not the signal 
itself, even right at the frequencies of maximum error. 

Conventional filtering, moreover, involves blocking 
frequencies in some range as strongly as possible, while 
passing others unaffected, and does so over a broad range 
of roughly similar installations. The method of the pre- 
sent invention, however, requires and involves the genera- 
tion of calibration data that reflects the detailed quan- 
titive characteristics of the supported mass and asso- 
ciated particular mechanical system of the specific 
installation . 

A further striking difference underlying the inven- 
tion is that it requires the output or outputs to have 
generally differing "filter" characteristics with 
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respect to each of the multiple inputs, determined by the 
specific mechanical system characteristics , Only in quite 
rare instances might two be the same, as previously 
pointed out . 

The novelty of the calibration procedure itself also 
appears to reside in several universal aspects. First, 
there is a large number of numbers to be determined — 
seventy two in the embodiment described — and they must 
be determined all at once, requiring at least the equiva- 
lent of a matrix inverse. A simple matrix inverse, how- 
ever, is generally impossible to use since, in practice, 
there is never enough of the right kind of data to define 
a unique and tailor-made calibration for a specific 
installation. One of many possibilities must thus, under 
the invention, be found which can be expected to perform 
satisfactorily* Finally, the calibration "input" con- 
sists, in accordance with the invention, of approximately 
specified disturbances, with such effected while the real 
system inputs with respect to what is to be measured, are 
allowed to remain at zero. 



Further modifications will also occur to those 
skilled in the art within the spirit and scope of the 
invention as defined in the appended claims* In addition 
to application of the several features of the invention 
independently or collectively to avoid the effects of 
inertial interference, these features, as before 
explained, may be applied to other types of force measure- 
ment systems, as well, including, but not limited to, 
weighing systems and the like. 



